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1. Introduction 

One dimensional (ID) quantum scattering theory is usually formulated for potentials 
that vanish asymptotically both for large positive and negative values of the coordinate 
x. It is well known that the degeneracy of the energy makes the full-line scattering 
problem somewhat more involved than partial-wave scattering on the half-line. 
Additional complications arise for step-like potentials, namely, when the potential tends 
to different constant values on both sides, 

lim V(x) = , lim V(x) = V > . 

x— oo x— >oo 

These conditions apply for example to electron collisions between different metals, in 
models of time-of-arrival measurement 0, 0, or in experiments with evanescent waves. 
In some cases it is enough to solve the Schrodinger equation numerically, subject to 
scattering boundary conditions, in order to obtain the transmission and reflection 
amplitudes. There are however applications where a formal theory of scattering 
is needed. By "formal theory" we mean the network of operators (Moller, S, T, 
and resolvents), which, together with their generic properties and relations, are used 
to describe the collision. These applications include the obtention of approximate 
analytical formulae, perturbative analysis, inverse scattering methods based on "two- 
potential formulae", kinetic theory, or the study of characteristic times ||. The work 
on the scattering theory of step-like potentials has concentrated on the inverse problem 
|TJ. |J, characterizations of scattering data for classes of potentials ||, zero energy 
limits 0], Levinson's theorem 0, and compact formulae for the evolution of states 
with initial support on one half-line || ||. This paper complements those mentioned 
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by focusing on the formal setting of the theory. In particular, we stress the fact 
that several partitionings of the Hamiltonian are possible, and work out, compare and 
relate the Moller operators and the corresponding Lippmann-Schwinger (LS) equations 
derived from them. Compact expressions of the asymptotic transmission and reflection 
amplitudes are given in terms of different potential-dependent matrix elements. The 
formalism is presented with "physicist's rigor". Its validity is in any case easily checked 
for cut-off potentials that deviate from the two asymptotic values and Vq only in a 
finite domain, [a, b], which is the case explicitly considered throughout. It is expected 
though that it will apply for other potentials as well, having smooth but sufficiently 
rapid decay. 

For completeness, and in order to introduce the relevant concepts and notation, in 
section II we present a lightnight review of Moller operators and Lippmann-Schwinger 
equations for potentials that vanish on both sides (the "ordinary case" hereafter), while 
some properties of scattering states of the Hamiltonian H for step-like potentials are to 
be found in section III. We discuss several partitionings of the Hamiltonian, together 
with the corresponding Moller operators and Lippmann-Schwinger equations in the 
following sections IV, V and VI. So as best to illustrate the differences among the 
formalisms we address the issue of the existence of Born's approximation in section VII. 



2. Moller operators for potentials that vanish on both sides 

In ordinary ID scattering the Moller operators Q±, defined by the strong limits 

Ot = lim e iHt/h e -iH t/n (1) 

link the actual state ip with its asymptotic free-motion reference states, (f>i n and <p ut, 
lim \\ip{t)- <pin/out{t)\ \ = , 

t—*— /+oo 

The operator Q + (respectively f2_) provides the scattering state by acting on the 
incoming (resp. outgoing) asymptote, (p in (resp. </w) , 

fi + /-|&n/aut(t)> = \lf>(t)), (2) 

for all t. 

The total Hamiltonian, H = Hq + V, is composed by a free motion Hamiltonian, 
Hq = p 2 /2m, that governs the motion of the asymptotes, and a potential operator, V, 
with a local coordinate representation = 8{x — x')V(x). The potential function 

V{x) vanishes as \x\ — ► oo, in such a way that the Moller operators in (|l|) exist. For 
concreteness, we shall in fact assume that V(x) vanishes outside the finite interval [a, b], 
with a < and b > 0. 

The infinite time limits in the definition of Q±, (|l|), may also be expressed with the 
alternative forms 

Ot= lim T e r°° dt'e^'e^'^e-^'l*. 

e^0±^ Jo 
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Inserting a resolution of the identity in momenta between Q± and \4> in (t)) or \4> ou t(t)), 
and integrating over t', there results 

/oo 
d pe -^A|p+/-)(p|0 in/out (O)), (3) 
-oo 

where we have introduced the (improper) eigenstates of H, with eigenvalue E p , 

|p±) = ft(£ p ± <0)|p> = \p) + —i-^V-jp), £ p = p 2 /2m. (4) 

± zO — if 

The states |p) (or \q), to be used in the following) are the usual plane wave states, 
(x\p) = exp(ipx/h)/v2nh (resp. (x\q) = exp (iqx /K) / \JlixK) . Vt{z) is a parameterized 
Moller operator (to be distinguished from the abstract ones in (|])) which, unlike Q±, 
can be applied to plane waves, and can be defined through 

Q(z) = 1 + G (z)f(z). 

In this equation, 

f(z) = V + VG(z)V (5) 

is the parameterized "T-operator" , or transition T operator, and the operators 

G(z) = (z-H)'\ and G (z) = (z - Hq)' 1 

are the resolvents for the Hamiltonians H and H respectively. Equation (|5|) is called 
the operator Lippmann-Schwinger equation. Expressions equivalent to (^) are obtained 
by using the operator Lippmann-Schwinger equation (]5p and the operator identity 

G (z)f{z)=G(z)V, 

which lead to 

^ = ^ + T-i I -r\ n f(E p ±iO)\p) = \p) + 1 ~ Vlp*). (6) 

hp ± l\) — Hq h p ± l\) — Hq 

Equations (f|) and (|]) are different alternative forms of the Lippmann-Schwinger integral 
equation for the states [p*). Note the structure of these states, composed by a free plane 
wave (incoming for \p + ), outgoing for \p~)) and a scattering part. The forms © are 
useful to determine the asymptotic behaviour of the states at large distances (for cut-off 
potentials this means x < a, x > b) since the matrix elements of Gq(E p ± iO) (the 
Green's function) are known, 

( r \ 1 | r '\ _ im r i(2mz)V 2 \x-x>\/h 

{l z-H 1 ' h(2mz)V* e 
In this expression the square root is defined with a branch cut along the positive axis. 
Using delta-function normalization (i.e., (p ± \p' ± ) = S(p — p')), the states behave outside 
[a, b] as 

( x \ sign( P )\ 1 y { ex P(W^) + r \p) exp(-ipx/h), x<a, 
h 1 / 2 \T l (p)exp(ipx/h), x > b, 
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. , . 1 | T r (— p) expUpx/fi), x < a. 

(x\ P -^) = —x\ V J ' \ . . _ '(8) 

n 1 I exp{tpx/n) + K {—p) exp{—tpx/n), x > b. 

Both in (|7|) and (|) p is a label for the energy. Let us first interpret the states in (|7|): 
for p > 0, there is an incident plane wave from the left, with wavenumber p/h, and 
R (p) and T l (p) are the corresponding reflection and transmission amplitudes for left 
incidence; on the other hand, if p < 0, there is an outgoing plane wave towards the left, 
with wavenumber \p\/fr, and the corresponding amplitudes are not properly related to 
"transmission" and "reflection" . However, since they are analytical continuations of the 
amplitudes for p > 0, the same notation is mantained. Similar considerations apply to 
the set of states described by (§). 

The particular form of the amplitudes T l,r (p) and R l,r (p) for potentials composed 
by square barriers is easily obtained by matching the wave function and its derivative 
at the edges. However, this procedure is useless in more general cases. Expressions of 
the amplitudes for the general case are obtained by comparing (0) and fl8|) with the 
coordinate representation of (P). In this way they can be related to on-the-energy-shell 
elements of the transition operators T(E P ± iO). We shall work out one case in detail, 
as a reference for later results. Assume p > and x > b. In 

dx'{x\G (E p + iO)\x')(x'\T(E p + i0)\p), 

-oo 

we can substitute \x — x'\ in the Green's function by x — x', since the support of 
{x'\T(E p +iO) \p) is necessarily restricted to be between a and b because of the dependence 
of T on V, see ([5]), and the finite support of V{x). Therefore, 

97TTO7 pip x /? 1 rOO ^ 

(x\p + ) = (x\p) - / dx'e- tpx/h (x'\T{E p + i0)\p) 



where 



h p 
(x\p) - ^^(x\p)T+ p , 



T^ = ( P \f(E p ±iO)\p'), \p\ = \p'\. 

Straightforward comparison with (0) leads to an explicit expression for T l (p). The rest 
of the amplitudes can be worked out similarly to obtain the following table 

np) =i-^^;r (p) , 

R\p) = _^TX P (P) , (9) 
p 

Note that time reversal invariance implies = T^ p , and therefore T r (p) = T l (p) = 
Tip). 
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3. Scattering eigenstates of the Hamiltonian for step-like potentials 

In the case of step-like potentials, the potential function V(x) does not go to zero both 
for positive and negative x, when \x\ —>■ oo. We shall assume in what follows that V(x) 
does indeed tend to zero as x — > — oo, and to Vo when x — > +00. In other words, we 
shall assume that V(x) equals Vg(x) = Vq6(x) plus some localized additional potential 
of finite support or that tends to zero sufficiently fast when \x\ — > 00. In such a case, the 
scattering part of the energy spectrum is doubly degenerate above Vo, as corresponds 
physically to incidence from one side or the other. Below Vo, however, there is only 
one linearly independent solution with an evanescent wave at x > 0. There may be 
bound states too, with energy Ej < 0. The resolution of the identity may be written in 
different ways, in particular as || 

l=Y,\E J )(E j \+ dp\p ± )(p ± \+ dp\p ± )(p ± \± dp\p ± ){p ± l 

■ J — oo Jpo Jo 

where po = (2mV"o) 1//2 and the states |p ± ), to be defined below, have as in the ordinary 
case an energy E p = p 2 /{2m). As pointed out above, p is a label of the energy. It can 
be positive or negative because of the degeneracy in energy. 

The states \p + ), with p > 0, have an incident plane wave of wavenumber p/h, and 
the states \p~), p < 0, a corresponding outgoing one, 

( x \jf&<p)\ 1 x ( exp(ipx/h) + R\p) exp(-ipx/h), x<a 
h 1 ' 2 1 T l (p)exp(iqx/h) 1 x > 

where q = (p 2 — 2mVo) 1 ^ 2 , with a branch cut that joins the branch points p = ±p , 
going slightly below lm(p) = 0. In this way the sign of q is the same as the sign of p for 
p 2 > Pq, whereas it becomes positive imaginary for — p < p < p . 

The states \p + ) for p < —po are defined by an incident plane wave from the right 
with wavelength —h/q(> 0), and states \p~) with p > p Q by an outgoing plane wave 
with wavelength h/q, 

{x \ p -^ ip)) = _J_ M V2 x jn-p)exp(zWfr), * < a 

h 1 / 2 \q J }exp(iqx/h) + R r (—p)exp(—iqx/h), x>b 

(always for \p\ > \po\). The factor (p/q) 1 ^ 2 is necessary for the proper delta 
normalization, that is, (p + \p' + ) = 5(p—p'), and the corresponding expression for the \p~) 
scattering states. As in the ordinary case, the arguments of transmission or reflection 
amplitudes are always positive for states \p + ), and negative for states \p~) independently 
of the sign of p. 

The S matrix elements are defined as the coefficients multiplying the outgoing plane 
waves when the incident plane wave is normalized to unit flux. When both channels are 
open (p > p ), the S matrix reads 

M- 1/2 < 

S(p) 




T\p) R\p) 
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One may also obtain these matrix elements from (p'~\p + ) by factoring out a delta 
function in the scattering energy. The unitarity of the S matrix, SS^ = 1, implies 
relations among the amplitudes, 

V(p)i 2 + i^ r (p)i 2 = i, 



q 

\R l (p)\ 2 + l\T l (p)\ = l, (12) 
V 



1\ T \ 
V 

'-T r ( P )R l ( P y + R r ( P )T l ( P y = 0. 



1 

For < p < po only one channel is open, the S matrix reduces to a number, R l (p), and 
unitarity implies 

r\p)r\pT = i. (is) 

All these equations, the set (12) and (13), are also valid for negative label p, thus 
providing relations for the amplitudes associated with \p~) states. 



4. Step-like potentials. Multichannel formalism. 

The straightforward application of the Moller operators of section II, based on the 
partitioning H = H + V, to step-like potentials is justified physically only for certain 
states. The key point is that H by itself only governs the asymptotic states that 
enter from the left (with incident positive momentum), or escape to the left (with 
negative outgoing momentum). So the Lippmann-Schwinger equations presented in the 
previous section (that is, eqns. (f|) and @), will only be applicable for {|p slgn ^^)}. 
It will prove useful to rename H as Hi = H , since it is the Hamiltonian associated 
with the "left" asymptotic channel. Correspondingly we define V\ = V, so that the 
total Hamiltonian is partitioned as H = Hi + Vi, and f2' ± = fi±. Similarly, the states 
|p-sign(p)^ > " s t ar t" (f or p < — p ) or "end up" (p > po) in the right, where 

the asymptotic Hamiltonian is H r = H + V$. We thus define V r = V — Vo, so that 
H = H r + V r , and the corresponding Moller operators 
Vi! = lim e iHt/h e -iH r t/n 

The asymptotic Hamiltonians have their own resolvents, 

G a (z) = 

z- H a 

where a = r, I is the subscript to indicate the channel. Notice that Gi(z) = Gq(z), 
using the notation of section II, whereas G r (z) = G (z — V ). Using the abstract 
Moller operators one may define parameterized ones, the corresponding LS equations 
thus taking the form 

| p si g n (p)) = lp) + £ l[Ep + S ign(^0]^b sisn ^> = 

= \p) + G[E P + sign(p)zO]V^>, (14) 
\p-^(p)) = \ qN ) + a r [E p - sign(^0]Kb- sign ^> = 

= \q N ) + G[E P - sign(p)iQ]V r \q N ), \p\ > \po\, (15) 
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where (x\qj^) = (p/hq) 1 / 2 exp(ixq/h). A noticeable difference with the ordinary case is 
that now the potential functions V a (x) are not localized {Vi(x) and V r (x) do not vanish 
for x > b and x < a respectively), so that the simple manipulations leading, for example, 
to (H), are not valid any more to obtain expressions for T l and T r . We cannot separate 
the exponential e l p\ x ~ x '\/ n [ n ^ x anc i x > dependent exponentials, and extract right away 
the x dependence. The separation can be done however to obtain R l and R r , which take 
the form 

R\p) = Z ^(-P\v l \p sien{p) ), (16) 
p 

R r (p) = ^^{qN\V r \ -p^), (17) 
p 

To obtain expressions for the transition amplitudes we rewrite the LS equations in terms 
of the potential of the other channel, see Appendix A, 

|p sign(p) ) =G r {E p + sign{p)0)V r \p sien{p) ), (18) 
| p - S ign( P ) ) = Q i{Ep _ sigri (p) )^b- sign ^) \p\ > \p \. (19) 

Since the potentials in fll8| ) and ( |i"9f ) vanish in regions of space different from the ones in 



(|l4l) and (|T^), we may now find the missing expressions for the transmission amplitudes, 

T\p) = ^^l(q\V r \p^), (20) 
Q 

T r (p) = ( -pM -p^). (21) 

P 

5. Jaworski-Wardlaw Moller operators 

In their study of the time spent by a quantum particle in a given spatial region ||, 
Jaworski and Wardlaw introduced two different asymptotic Hamiltonians for incoming 
and outgoing asymptotes, 

H m + V F. 

2m 

— - p 2 

H ou t = h VqF+ 

2m 

where F_ and F + are complementary projectors, F_ + F + = 1, over negative and positive 
momenta respectively, 

^ r±oo 

F ± = ± dp\p){p\. 
Jo 

Correspondingly, they defined Moller operators 

Q J JY_ = lim jHt/h e -iH in/out /h_ / 22 x 

Note that, as in the previous section, two different partitionings of the Hamiltonians 
are required, one for each Moller operator. They are however not based on right/left 
channels, but on a distinction between incoming and outgoing states. The physical 
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reason for these definitions is clear: the positive momentum part of the incoming 
asymptotes travels on the lower level at long negative times, whereas the negative 
momentum parts travels on the upper level. The outgoing asymptotes behave in the 
opposite way, with positive momenta on the upper level and negative momenta on the 
lower level at large positive times. 

We shall now extend this formalism to produce the asociated Lippmann-Schwinger 
equations. First it is convenient to introduce a delta-normalized eigenbasis for H in and 
H out , (explicitly, (in(p)\in(p')) = S(p — p'), and similarly for \out{p))) 

( jpx/h p>0 

\\p/Q\ 1/2 ^ iqx/n P<-Po 

{ Ip/ql^e^ p>p 

\ e ipx/h p<0 
so that 



(x\in(p)) = h x ' 2 x 
(x\out(p)) = h~ x l 2 x 



Hi n / ou t\in/out{p)) = E p \in/out(p)). 
Aside from the ordinary (momentum) resolution of the identity, 1 = J^dp \p)(p\, 

1= / dp \in(p)) (in(p)\ + / dp \in(p))(in(p)\ = 

J-oo Jo 
rO roo 

dp\out(p))(out(p)\ + / dp \out(p)){out(p)\. 

The connection between the abstract Moller operators (^) and Lippmann-Schwinger 
equations for eigenstates of H follows now closely the steps from (fj) to @, but making 
use of the above resolutions of the identity. We thus find 

dp\p + )(tn(p)\<f) m (t))+ / dp\p + )(tn(p)\<p m (t)), 

-oo JO 
rO roo 

\4>(t))= dp\p-)(out(p)\<j) out (t)} + dp\p-)(out(p)\(j) out (t)), 



with 



and 



\p+) = \in(p)) + G in (E p + iO)f in {E p + iO)\in{p)), (23) 

\p~) = \out(p)) + G out {E p - tO)f out {E p - iO)\out(p)), (24) 

Gi n / ou t(z) = (z — Hi n / ou t) , (25) 

Ti n /out( z ) — Vin/out + Vin/out G(z) Vinjout > (26) 



corresponding to the two partitionings of the Hamiltonian, H = H in + V in = H out + V out , 
where 

VL = V-V Q F-, and V out = V - V F + . 

However, the potentials Vi n and V ou t are not localized. They do not vanish as x — > oo 
and this leads to similar problems to the ones encountered before when searching 
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for expressions for the transmission amplitudes. They are actually more severe now 
because these potentials do not have a semibounded support; in addition, the zeroth 
order Green's functions, which can be explicitly obtained by integration in the complex 
momentum plane, are cumbersome to work with, 

(x\G in (z)\x') =(x\— ^+\ x ') + ( x \ ^—^\ x >), (27) 



z — H z — Vq — H i 







(x\G out {z)\x>) = (x\ ^—~W) + (x\-^\x'). (2f 

z-Vo-H z- H 



The summands are particular cases of 



F « \ x ') = A£sign(x - x') + 9[£(x - x')}(x\Go(()\x'), 



(-Hi 



with £ = ± and 



y = (2m() 1/2 \x-x'\/h, (29) 

where the square root with positive imaginary part is taken. 

The scattering states defined through ( p3|) and ( pl|) are the same as those defined 
previously by the LS equations of the multichannel method presented in the previous 
section. In order to check the veracity of this statement, it is convenient to use the 
identity Gi n / ou t(z)T in / ou t(z) = G(z)Vi n / out , together with the forms of the resolvents G 
given in ([y]) and (|I5|) , and the defining expressions of the different potential operators 
involved. 



6. Pure-step Hamiltonian as zeroth order 

In this section we shall study one more possible partitioning of the Hamiltonian, by 
considering the Hamiltonian H s = H + V e for the pure step potential Vq = V 8(x), as 
the zeroth order term for the complete Hamiltonian, 

H = H S + V S . 

In other words, the total potential energy is decomposed into the pure step potential part 
and a localized part, as V = Ve + V s . It is easy to compute two different eigenbases of H s , 
whose elements are \pf) respectively (labeled by p as before). Their explicit expression 
lends itself to identification of transmission and reflection amplitudes by comparison 
with expressions (pTOj ) and (|Tl|): 

q+p q+p 

TZ(p) = K( P ) = ^ (30) 

p + q p + q 
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Green's function for H s is also known exactly [JTO 



r e ±t|p||a;-a:'|/7i _|_ r±e =F*|p|(aH-a;')/'»l 



(x\G s (E p ±iO)\x') = ± 



m 
ih 



( 1 

¥\ 

1 

¥\ 

J_ . e ±i(fj,±x'-\p\x)/h 



t±e 



±i(/i±x— \p\x')/h 



\P\ 

_}_\ e ±iH±\x-x'\/h _ r ^ e ±i[t±(x+x')/h-i 



x' < 0, x < 
x' < 0, x > 
x' > 0, x < 
x' > 0, x > 



(31) 



where -E p = |p| 2 /2m, 



2\p\ 



\P\ +M±' 



and 



11/2 



[2m(E p - V Q 
±t[2m(V -E p )] 1 / 2 , 



E p > V 

E„ < Va 



A first advantage of this decomposition is that the state is governed asymptotically by 
H s both before and after the collision, to the right and to the left, so that the physically 
meaningful Moller operators can be defined, as in the ordinary case, by the two limits 
of a unique operator expression, 

Q s ± = lim e iHt/h e -iH.t/h^ 



t— >=poo 



which amounts to a formal simplification with respect to the partitionings of the two 
previous sections, and absence of extra indices. Analogous steps to those leading to @, 
with the decomposition of unity in the basis of H s , provide us with 

~P0 fOO 

dp\p + )(pf\<p m (t))+ / dp\p + )(pt\(t>in{t)), 



where 



and 



\m) = 

\m) = I dp\p-)(j)J\(f> tmt (t))+ I dp\p-)(p- s \<p out (t)) 

J — oo J po 

|p ± ) = \pf) + G S (E P ± i0)f s (E p ± i0)\ P f) 



(32) 



T s {z) = V s + V s G{z)V s . 

A second advantage of this decomposition is that the potential function V s (x) = 
V(x) — ^^(x) is now localized. One may thus obtain easily the explicit x-dependence of 
the coordinate representation of (B2|) and identify expressions for the transmission and 
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reflection amplitudes in terms of the localized potential, 

T\p) =Tj(p)-^p(j) ' (p; si ^)|f s [^ + sign(p)^0]b^), 

T r ( P ) = T:(p) -^(^Y (- P ;^\f s [E p + sign(p)iO]| , 
R\p) = R l s (p) - ^(- p7 ^(p)\f s [E p + sign(p)^0]b s ;^)) (33) 

Aj 7T 7 77 ? 

iT(p) = i£(p) - ——(j, 7 ^(p)\T s [E p + sign(p)iO]| - p^W) , 

The time-reversal antiunitary operator G changes the sign of Q±, 0Q± = Q T Q, as in 
the ordinary case. From the time reversal invariance of the Hamiltonian it follows that 
(p\T s [E p + sign(p)iO]\p') = ( — p\T s [E p + sign(p)zO]| — p') (on the energy shell), so that 
the transmission amplitudes are related by T r (p) = {q/p)T l (p). 

The agreement with the previous compact expressions ( p!6| - |T7|) , and ( f20| - |2l| ) is found 
by using ( p^ - |T7|) and ( ^0| - |2T| ) themselves for the step potential Vg, and the following non 



trivial generalizations of the standard "two-potential" formula to the two partitionings 
of the multichannel formalism (see Appendix B), 

(p\V\± j9- sign(p) ) =(p\V e \± p; si ^) + (pf gn{p) \V s \ ± p- si ^), 

(q N \(V - V )\ ±p si ^) = (q N \(V e - V )\ ±^ ign(p) ) + ( P ; sign(rt \V S \ ±p sign M>. 

The use of the bra-ket notation, while standard and very convenient most of the 
time, requires some greater attention than usual to describe adequately the evanescent 
case, when q = ry, 7 > 0. Irrespective of the value of p, (q\x) should always be 
interpreted as h~ x l 2 exp(— iqx). Similarly, (pf\x) should first be written for real q and 
then continued analytically. 

7. Born approximations 

As an example to illustrate the differences of the three described formalisms we shall 
obtain the Born approximation of the reflectance \R l (p)\ 2 for the potential, 

V = V e + V 1 6(x). (34) 

The exact result, 

, = p-q- (2imV 1 /h) 
{P) p + q+ (2imVi/h) ' 

may be obtained using (pOj) , (|33|), 5(x) = |0)(0|, and 

f ( ) = Y±BE 

s{z) 1-^(01^)10)' 

or alternatively by straightforward computation. 

We will now calculate the different Born approximations by retaining only the 
terms linear in the potential corresponding to each partitioning of the Hamiltonian. 
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To be more precise, we will look at the Lippmann-Schwinger equation for \p + ) in each 
approach, and retain terms of first order in the potential, that is, to first order in the 
difference between the total hamiltonian and the incoming asymptotic hamiltonian of 
reference. The resulting wavevector will be examined in the position representation for 
x < 0, and the result compared to flIU|) to extract R l {p). In fact this last step is not 
necessary for the multichannel (MM - sec. IV) and the localized potential (LP - sec. 
VI) approaches, since we have already carried out this comparison in an exact manner 
(see (0) and (|3~3f)). Notice that we have indeed checked that in all three approaches we 
obtain the same scattering state \p + ). 

In the multichannel method (MM), see (fL6[), the first order term in Vi (which in 
this case is Vi — Vq + V\5(x)) is 

/ x —2nmi. T ^ , m(V — 2ipV 1 /h) 

R l Born-MM(p) = ( ~ P Vg + V l 5{x)\ P ) = f (35) 

p 2p z 

The analysis to second order is more delicate, involving limits (as in G(E p + i0)), but it 
reveals that the singularity in p = for the MM formalism actually worsens (it becomes 
of the form p~ A ) . This was only to be expected, given the non locality of the perturbing 
potential in that case, which produces infrared singularities to all perturbation orders, 
which can only be resolved by a complete resummation of all terms. 

We could also examine R r (p), to first order in V r = Vg + V\8{x) — Vo, which in 
this case can be obtained from R l (p) by substituting p for q, and viceversa. This recipe 
actually holds for the Born approximation in the MM formalism, which leads to the 
result that the reflectance diverges for \p\ — ► p . 

A Born approximation in the "in/out" formalism of section V is much more 
problematic: for x < and p > 0, the first order in Vi n of (|23"D is 



x\p + ) {1) = (x\in(p)) + [ dx' (x\G in (E p + iO)\x')(x'\V in \in(p)). 



By substituting (x\Gi n (E p + i0)\x'), see (|27f) and (|29|) , and taking the limit x — > — oo to 
eliminate transient terms, 



(x\p + ) {1) = -j= 



e i P x/n + ( P-Q + mVA e - iqx/h 



2q ihq 

To this order, this approach provides a physically meaningless reflected wave with a 
momentum smaller than the incident one. This indicates that we do not recover in 
this manner a sensible approximation to the reflectance. Finally, the localized potential 
(LP) approach of the previous section gives, to first order in V s , 

RBorn-Lp(P) = — (~Ps 0) pj) = . 36 

p + q p {p + q) 



The results of (|35| ) and (^g) are compared in figure 1, which clearly demonstrate 
the computational advantage of the localized potential approach, which starts from 
much better adapted initial functions. In particular it is relevant to note that the 
LP approach detects the change of regime in the reflectance due to the energy falling 
below the asymptotic level, which the multichannel formalism cannot even suspect in a 
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Figure 1. Exact reflectance (solid line), first order Born approximations for the 
localized potential approach (dots) and multichannel method (dashed line), and second 
order approximation for the multichannel method (squares) versus p. The potential is 
given in (|34|). Vq = 1 a.u., V\ = 0.01 a.u., m = 1 a.u. 

perturbative scheme. In other words, the fact that only one channel is open, and (pf) 
must hold is overlooked by the perturbative expansion in the MM scheme, while there 
is a sharp change in behaviour of the perturbative expansion in the LP scheme from the 
one channel to the two channel case (even though ( |T3"D does not generically hold if we 
restrict ourselves to a finite number of terms). 

8. Discussion 

Given the simplicity of one dimensional step-like potentials, we could not fail to provide a 
complete formal scattering theory for them. However, in pursuing this objective, we have 
met several interesting novel aspects with respect to ordinary scattering. Among them, 
the existence of different, all somehow "natural", partitionings of the Hamiltonian is an 
important one, since it leads to different formal frameworks. Working out the details is at 
the very least laborious, frequently tedious, and we hope that our compact presentation 
and focus on the final results may save some time and help the practitioners to avoid 
pitfalls. With respect to the three possible methods described, the in/out-formalism 
has some elegance, and this was historically our first choice. However, the zeroth order 
Hamiltonians, non localized potentials, or Green's functions are not easy to deal with. 
This lead us to look for other possibilities. Certain manipulations may benefit from 
the condensed forms of transmission and reflection amplitudes obtained following the 
multichannel method, but in a practical calculation, the localized potential approach 
will be generally preferable. It is also the most economical presentation since it reduces 



Step-like potentials 14 
in half the number of equations needed, and is also the closest to ordinary scattering. 
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Appendix A. Alternative forms of Lippmann-Schwinger equations 

We give an example of the obtention of the alternative LS equations in ( |j"8|) and fll9|) . 
Using 

Gi(z) = G r (z)[l-V G l (z)] 

and V r = Vi — V , ([14]) for p > may be written as 

\p+) = \ p ) + G r (E p + iO)V r \p + ) + V Q Gt(E p + zO)[l - G r (E p + iO)V r ]\p + ). 

Acting with the operator in parenthesis on \p + ), using G(z) = Gi(z) +Gi(z)ViG(z), and 
(P^), one finds (when operating +i0 must be kept as a small imaginary number) 

\p+) = \ p ) + G r (E p + iO)V r \p + ) + V G r (E p + iO)\p), 

but the third term cancels the first one by acting with G r on \p), so that fll8f) (for p > 0) 
is obtained. One may proceed similarly for the other cases. 

Appendix B. Two potential formulae 

In this appendix we shall obtain one of the two potential formulae used in section VI. 
The other cases may be obtained similarly. Assume that p < 0. Then, 

(p\V\±p + ) = (p\V e + V s \±p + ) 

= (p\(V e + V s )[\ ±pf) + G S (E P + iO)T s (E p + z0)| ±pf)] 

= (p\V e \ ±pt) + (p\V e G s (E p + tO)f s (E p + i0)\ ±p+) + (p\V s \ ± 

= (p\Ve\±pt) + (p;\V s \±p + ), 

where we have used ( |32"D and (p!4]), the last one particularized for the pure step potential. 
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